Magnetism, Spin-Orbit Coupling, and Superconducting Pairing in UGe2 
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A consistent picture on tlie mean-field level of the magnetic properties and electronic structure 
of the superconducting itinerant ferromagnet UGe2 is shown to require inclusion of correlation 
effects beyond the local density approximation (LDA). The "LDA+U" approach reproduces both 
the magnitude of the observed moment, composed of strongly opposing spin and orbital parts, and 
the magnetocrystalline anisotropy. The largest Fermi surface sheet is comprised primarily of spin 
majority states with orbital projection rai=Q, suggesting a much simpler picture of the pairing than 
is possible for general strong spin-orbit coupled materials. This occurrence, and the quasi-two- 
dimensional geometry of the Fermi surface, support the likelihood of magnetically mediated p-wave 
triplet pairing. 



The coexistence of superconductivity with magnetism 
has recently re-emerged as a central topic in condensed 
matter physics, due to competition between ordering of 
magnetic ions and superconducting pairing in borocar- 
bides to unprecedented coexistence of magnetic or- 
der with substantial ferromagnetic component with high- 
temperature conductivity [0-H and to observation of 
magnetic superconducting transitions near quantum 
critical points in f electron systems such CeRh2Si2 
Celna |], CePd2(Si,Ge)2 §, and CeRhlns g. Super- 
conductivity below IK in the limited pressure range (1 
to 1.6 GPa) was recently observed coexisting with strong 
ferromagnctism in a high purity single crystal UGe2 [p|JlO| 
adding yet another dimension to this unanticipated phe- 
nomenon of phase coexistence. This superconducting 
phase is found within the ferromagnetic phase and dis- 
appears in the paramagnetic region, strongly suggesting 
the pairing mechanism is magnetic in origin. 

There are no other known examples of superconductiv- 
ity coexisting with strong ferromagnctism (average mo- 
ment comparable to Ni, say). This phenomenon neces- 
sarily requires triplet spin pairing (or more precisely, an 
odd symmetry real space pair) since the exchange split 
bands (and the very different Fermi surfaces) preclude 
singlet pairing. The normal state symmetry itself is low 
in spite of the relatively simple crystal structure of UGe2. 
Magnetism breaks time reversal invariance, leaving the 
normal state symmetry consisting of the product of the 
magnetic space group and gauge invariance [Q x C/(l)], 
the latter of which is always broken by the onset of super- 
conductivity. The possible order parameter symmetries 
are very limited. 

Recent experiments on single crystals [ pdp^ ] indicate 
UGe2 (previously thought to crystallize in the Cmcm 
ZrSi2 type structure) actually has the base-centered or- 
thorhombic ZrGa2 crystal structure {Cmmm). The 
structure, shown in Fig. 1, can be viewed as consisting of 
antiphase zigzag chains of U atoms running along the a 
direction and lying within the a — b plane; however, inter- 
chain and intrachain separations are comparable. Each 
U is tenfold coordinated by Ge. Importantly, the struc- 
ture possesses inversion symmetry. Single crystal magne- 
tization measurements |0] and neutron powder diffrac- 



tion measurements [ |l2[ both yield a coUinear magnetic 
structure with the ferromagnetically ordered magnetic 
moment of 1.42 fiB- (We quote moments per formula 
unit, i.e. per U atom.) The Curie temperature Tc — 52 
K at ambient pressure decreases under pressure, vanish- 
ing at 1.6 GPa. Around 1 GPa, Saxena et al. ||l^ have 
found that UGe2 becomes superconducting while remain- 
ing strongly ferromagnetic {M « 1 /xb/U) providing 
a novel example of coexistence of the superconductivity 
with strong ferromagnetism. 

Magnetic measurements [^,0 yield very strong 
magneto-crystalline anisotropy in UGe2 with easy mag- 
netization axis along a (the shortest crystallographic axis. 
Fig. 1). It was found jl^ that even in high magnetic field 
(up to 35 T) it is impossible to saturate magnetization 
along the hard b and c axes. Such highly anisotropic 
behavior is typical of magnetic 5/ electron materials. 
Susceptibility measurements pO[ | above Tc yield an ef- 
fective paramagnetic moment me// ~ 2.7 /is/U which 
differs from atomic value of 3.62 /j.^ for atomic con- 
figuration, indicating substantial 5f hybridization with 
conduction states. Moreover, the carriers with the high 
cyclotron masses of 15-25 mo were observed in dHvA 
experiments 15 1 suggesting itinerant but strongly corre- 
lated 5f-electron states. UGe2 is however not a real heavy 
fermion system since its electronic specific heat coefficient 
C (T) /T = 7^35 mJ /K^mol is about ten times smaller 
than in conventional heavy- fermion U-compounds jl^ . 

The electronic and magnetic structure of UGe2 is a 
rather open issue. Yamagami et al. p6[ | reported the re- 
sults of relativistic augmented-plane-wave calculations, 
but without interplay between magnetism and spin-orbit 
coupling (SOC) and using the old crystal structure. The 
scope of this paper is to provide a microscopic picture 
of electronic structure and magnetism of ferromagnetic 
UGe2 upon which models of the superconductivity can 
be based. This picture leads to specific indications of the 
pairing character in UGe2. 

We first applied conventional band theoretical meth- 
ods in the local density approximation using a relativistic 
full-potential linearized augmented plane wave method 
I prt . The results did not reproduce the observed ground 
state moment ||lj, |l|]. The calculated spin magnetic mo- 



1 



ment {Ms = 1.57 /is) is canceled by the orbital moment 
arising from spin-orbit coupling {Ml = -1.84 |18|) 
yielding the net moment of \M\ = 0.27 /is, five times 
smaller than the experimental value of 1.4 /i^- The cause 
is due to the oversimplified treatment of correlation ef- 
fects, as is often seen in the application of LDA to f elec- 
tron materials. 

We therefore account for the on-site atomic-like corre- 
lation effects beyond LDA by using the LDA-I-U approach 
in a rotationally invariant, full potential implemen- 
tation. Minimizing the LDA+U total energy functional 
with SOC treated self-consistently [|l^ generates not only 
the ground state energy and spin densities, but also effec- 
tive one-electron states and energies (the band structure) 
that provides the orbital contribution to the moment and 
Fermi surfaces. The basic difference of LDA-I-U calcula- 
tions from the LDA is its explicit dependence on the on- 
site spin- and orbitally resolved occupation matrices [ pO[ . 
Since the LDA-I-U method is rarely applied to metals, the 
appropriate values of the on-site f electron repulsion U 
and exchange J constants are not known. Our values 
of U — 0.7 eV, J — 0.44 eV were chosen to reproduce 
the ground state magnetic moment Mj = Mg + Mi. 
Mainly U was varied, but it seemed to be of some impor- 
tance to reduce J slightly from its atomic value of 0.55 
eV |2^ to keep U from becoming unreasonably small. 
The resulting total uranium f-electron occupation (2.8 e) 
is close to but clearly less than the configuration. The 
value of Coulomb U differs from the "atomic-like" value 
oi U — 2 eV which was derived for U compounds in 
Rcf. 1^. The reduction of the Coulomb U to 0.7 eV in 
order to obtain the correct value for the magnetic mo- 
ment reflects the partly itinerant character of the U 5f 
states and accounts for the metallic screening of the U 5f 
states in UGe2. 

To probe the magnetic ground state we calculated the 
energy of one antiferromagnetic configuration in addi- 
tion to ferromagnetic, where the two U atoms in the 
unit cell (in the "chain" in Fig. 1) had antialigned mo- 
ments. Nearest neighbors between chains there were also 
antiparallel. The ferromagnetic alignment is 32 meV/f.u. 
lower in energy, consistent with the observed ferromag- 
netism. 

To study the magnetic anisotropy we performed the 
total energy calculations for the magnetic moment along 
each of the a, 6, and c axes. The orthorhombic crystal 
symmetry was reduced in order to preserve the projection 
of the orbital moment along the spin direction. To obtain 
convergence for the total energy differences we used 275 
k-points in l/4th of the BZ for the magnetization along a 
and b, and 250 k-points for c orientation; these meshes are 
equivalent to 1000 k-points in the BZ. The total energy 
differences yield the magnetic anisotropy energy and are 
shown in Table I. The calculations produce easy magne- 
tization axis a in accord with the experiment. There is a 
pronounced anisotropy for the both spin M^ and orbital 
Mi magnetization reflecting the interplay between the 



strong spin-orbit coupling (the SOC constant ^ = 0.22 
eV) and the large exchange splitting. The variation of 
the spin and orbital moments with orientation is 7-10% 
(Table I), and manifests itself in strong, qualitative dif- 
ferences in the band structure and Fermi surfaces. Thus 
it is crucial to use the easy axis orientation, and hence- 
forward we discuss only this case. 

This strong magnetic anisotropy (2-3 orders of magni- 
tude larger than in 3d ferromagnets) reflects the strong 
spin-orbit coupling for 5f-states and is consistent with 
the reported experimental data for the MAE in the ura- 
nium compounds p^ . It explains qualitatively the dif- 
ference in the measured magnetization values for single 
crystals (1.42 (Xb) and polycrystals (1.07 ^b |@, 0.5 fj,B 
(|l4|): since the magnetic moment is kept by magnetic 
anisotropy along the easy magnetization direction (a) in 
each grain of the polycrystal, the average measured mo- 
ment should be reduced but should exceed one third of 
its single crystal value. 

The complex band structure contains 7 orbitals x 2 
U atoms X 2 spins = 28 primarily U 5f bands within a 
3 eV region spanning the Fermi energy Ep- They are of 
mixed spin character due to SOC and contain only minor 
Ge 4p character, which is mostly below and above the U 
5f bands. Bands crossing the Fermi level have predomi- 
nantly U 5f character and mostly spin majority character. 
These bands have low dispersion along ky {b direction) 
indicating an a — c plane quasi-two-dimensional (quasi- 
2D) character of the UGe2 electronic structure near Ep- 
The exchange splitting, to the extent that an average can 
be defined when there is strong SOC, is about 20 times 
as large (1-1.5 eV) as suggested by Yamagami et al. |l^] 
based on unpolarized calculations. 

Since the LDA-I-U potential is formulated in terms of 
TOs and mi rather than the total moment m^, it is nat- 
ural to perform a spin- and orbital- resolution {ms;m{\. 
The combination of strong spin polarization and large 
SOC result in a remarkably clean {ms]mi\ separation as 
is evident in Fig. 2. The peak in the vicinity of E^^^ is 
formed by {T;0} and {|; 1} states coupled by spin-orbit; 
N{EF)/i.Vi. = 1.7 and 0.7 states/eV respectively. The 
total N{EF)/tn. of 5.5 states/eV is mainly (75-80%) due 
to U 5f contributions. The measured electronic specific 
heat coefficient 7 = 35 mJ/K^ mol |14| corresponds to a 
dressed value N*{Ep) = 15 states/eV, indicating a dy- 
namic enhancement N*{Ep)/N{Ep)— 2.7 arising from 
magnetic fluctuations with possible contributions from 
phonons and charge fluctuations. 

Photoelectron (XPS-BIS) studies of UGca |2|] re- 
ported a pronounced broad peak in XPS spectrum « 1 
eV below the Fermi level. Our calculated U fDOS (Fig. 
2) has a double peak due to {T;w/ = —3,-2} located 
0.3-0.4 eV below the Fermi level. This difference may be 
due to neglect of dynamical effects in our calculations but 
may also reflect surface state contributions in the data. 
In the measured BIS spectrum, which is less surface sen- 
sitive, there is a kink around 0.5 eV and a main peak 
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at 1.3 eV above Ep. The LDA+U calculations repro- 
duce both features: the structure at 0.5 eV arises from 
{t; — 1} character and the much bigger peak at around 
1.5 eV is due to {t; +3} and {J,; 0, 1} states. 

The calculated Fermi surface is complex and consists 
of three sheets. The largest and most interesting sheet is 
shown in Fig. 3. The quasi-2D surface has weak disper- 
sion along ky, i.e. hopping between chains is small for 
this band. In spite of (but actually due to) the strong 
SOC, this surface has simple character: predominantly 
{t;0} with some mixture of {|; 1}. If, in a first ap- 
proximation, this small opposite spin character can be 
neglected, one can consider pure spin {t;0} states and 
hence simpler models for pairing. In this limit the prob- 
lem reduces to the one of so-called "single spin pairing" , 
which however has not been studied specifically for 
the low orthorhombic symmetry of UGe2 . 

It was suggested some time ago [|5| that supercon- 
ducting p-wave (triplet) pairing for equal spin states 
can appear due to the longitudinal magnetic fluctua- 
tions. The roughly nesting portions of quasi-2D sur- 
faces should promote strong magnetic interactions be- 
tween spin-majority carriers located periodically in the 
space along magnetization direction a. The rapid col- 
lapse of the moment between 1 GPa and 1.6 GPa |l^ is 
also consistent with longitudinal magnetic fluctuations. 
The distinctive electronic structure that we find suggests 
that a non-relativistic model can be used to describe the 
magnetically mediated superconductivity in UGe2. 

It must be recognized that the results presented here 
apply only to the zero pressure and temperature. At the 
peak of the superconducting phase (~ IGPa) the moment 
is reduced by one-third and the lattice constants and in- 
ternal coordinates will be changed. Consequently the ex- 
change splitting will decrease and band widths will be in- 
creased, altering the Fermi surface from the one we have 
presented. The general features of the {ms\mt) fDOS 
separation in Fig. 2 are expected to be robust. 

To summarize, accounting for some correlated behav- 
ior (beyond what is included in LDA) using the LDA-I-U 
approach in a very general implementation including 
SOC, we have provided a microscopic picture of the 
electronic and magnetic character of UGe2 that is es- 
sential for beginning to understand the coexistence of 
superconductivity with strong ferromagnetism. Our re- 
sults are in accord with several important experimental 
observations (magnitude and direction of the moment, 
BIS spectra), and suggest a relatively simple picture of 
the U 5f bands arises in spite of the 28 f bands that 
span Ep. The calculated Fermi surface has a quasi-2D 
sheet of predominantly majority spin character, support- 
ing conjectures of magnetically mediated triplet super- 
conductivity in UGe2. The coexistence of superconduc- 
tivity with magnetic order is quite distinct from that in 
RuSr2GdCu2 08, with its separate magnetic supercon- 
ducting layers, because in UGe2 the magnetism and su- 
perconductivity clearly seem to be due to the same U 5f 
electrons. 
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TABLE I. The U 5f state occupation (n/), spin (Ms), or- 
bital (Ml) and absolute value of net (|Mj|) magnetic mo- 
ments (/xb) for three directions of magnetization. The mag- 
netocrystalline anisotropy energy (MAE in mRy/f.u.) is the 
total energy difference MAE = Eij(c) ~ Ea for the magne- 
tization directed as indicated. 
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\Mj\ 


MAE 
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2.81 


1.52 


-2.98 


1.46 
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2.80 


1.41 


-2.85 


1.43 


0.67 
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2.81 


1.62 


-3.18 


1.56 


0.55 



Exp. |jjl2 



1.42 




Ge 



FIG. 1. The base centered orthorhombic Cmmm crystal 
structure of UGe2. The volume shown includes two primitive 
cells. 
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FIG. 2. Spin and orbitally resolved U 5f state density of 
states for easy axis a orientation: spin-| is shown in black full 
line, spin-| is in blue dashed line. Other orientations of the 
moment give distinctly different results. 



2.5 



4 




FIG. 3. The Fermi Surface for UGe2 that is discussed in the text. The zone center is at the far corner of the plotted regions, 
and the distances plotted correspond to ^ with axes directed as in Fig. 1. The large surface is primarily {nis =T; m; — 0} 

character, is quasi- two-dimensional and shows indication of a nesting feature around (0.452^^,0,0). 
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